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The transverse and longitudinal current correlation functions are evaluated in liquid and amor¬ 
phous sodium by computer simulation. The study of the corresponding spectra as a function of the 
wavevector k allows the evaluation of a dispersion curve for their peak position. The results are 
compared with recent experimental findings [PNAS 107 , 21985 (2010)], obtained by a new analysis 
of S{k,ij) measured by inelastic X-ray scattering in liquid and poly crystalline sodium. A substantial 
agreement between experimental and computer simulation results is found. The analysis of the line 
widths supports the Vibration-Transit theory picture of the dynamics in liquids. 
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The use of modern X-ray sources and spectrometers 
has allowed in recent years the measurement of the spec¬ 
trum of density fluctuations with surprising great accu¬ 
racy at very small frequencies, a range where neutron 
spectroscopy has no access. In particular the Brillouin 
frequency range could finally be explored. One could 
study in great detail the departure from the hydrody¬ 
namic regime, where the density fluctuations are sam¬ 
pled over a distance much larger than the interatomic 
spacing. In particular liquid metals have received great 
attention at the very early stage of the new experimen¬ 
tal technique: results have been thoroughly reviewed W 
Scopigno, Ruocco and Sette in an exhaustive paper [ij, 
where an extensive list of references to the subject can 
also be found. 

Parallel to the experimental effort there has been the 
attempt to present a comprehensive theoretical interpre¬ 
tation. This has been achieved by resorting to the Mori- 
Zwanzig ii approach. The method relies on the pro¬ 
jection operator formalism to arrive to a second order 
differential equation for the relevant correlation function. 
This correlation function contains the first order memory 
function, for which parameters have to be found by fitting 
experimental data. An interpretation in terms of under¬ 
lying physical processes is provided by recognizing the 
effects of short time (free particle) dynamics and struc¬ 
tural relaxation. It is useful to remind that in this con¬ 
text the adopted projection operator method turns out to 
be equivalent to the assumption that the hydrodynamics 
point of view can be extended at the microscopic level 
if transport coefficients (diffusion, viscosity etc.) are re¬ 
garded as both space and time dependent . In particu¬ 
lar the memory function approach never involves quanti¬ 
ties like generalized viscosities (shear and bulk), since the 
experimental spectra concern only the longitudinal cur¬ 
rent, whose spectrum is directly related with the experi¬ 


mentally available dynamical structure factor S{k,u!). In 
the hydrodynamics framework in fact longitudinal and 
transverse currents are orthogonal. 

An alternative approach, known as Vibration-Transit 
theory (V-T), has been more recently proposed by Wal¬ 
lace y. A liquid is treated as a disordered crystal with 
the same theoretical techniques valid for such a struc¬ 
ture. The theory not only gives a good account of most 
of the thermodynamic properties of liquid metals already 
in its zero order version, but it also provides a tool for 
describing the dynamical effects 0. In this framework 
the system is pictured as vibrating in what are referred 
to as random valleys, plus instantaneous transits among 
neighboring valleys. In this theory the longitudinal and 
transverse currents do not play a significant role even if 
their expressions can be calculated and indeed reported 
in detail 0- 

Molecular dynamics (MD) calculations have played an 
important role in confirming the physical interpretation 
since the early stages of computer simulation techniques. 
In a seminal paper Rahman showed that liquid rubid¬ 
ium supports propagating density fluctuations for wave 
number k < 1.2A~^ and found good agreement with neu¬ 
tron scattering results . Similarly the direct evaluation 
of the spectrum of the density fluctuations by computer 
simulation confirmed the accuracy of the V-T theoreti¬ 
cal results for liquid sodium 0, when compared with 
experimental data [Hj. 

The fact that density fluctuations, as measured by neu¬ 
tron or e.m. scattering, reflects only the longitudinal 
current dynamics, became deep-rooted to such an extent 
that in dealing with the study of transverse current corre¬ 
lation by computer simulation, it was often reported that 
such a quantity cannot be measured by any experimental 
method. However the serendipity of molecular dynamics 
techniques opened new routes for the most comprehen- 


2 




Frequency (THz) 




2SSssi«», 


Frequency (THz) 



Frequency (THz) Frequency (THz) 


FIG. 1. Spectra of the transverse current at different wavevec- FIG. 2. Same as Fig. [T]but for the longitudinal current, 

tors k, shown on each graph in A~^ units. 


sive interpretation of the physical processes underlining 
the scattering events. In studying the dynamical proper¬ 
ties of water, modeled by a new polarizable interaction 
potential, Balucani et al. 12| pointed out that the spec¬ 
trum of the longitudinal current shows the presence of 
two separate peaks at wavevector k > 1.0lA“^, one of 
them being coincident with the peak of the transverse 
current at similar wavevectors. Moreover, in the same 
range of wavevectors, the transverse current spectrum 
was found to display peaks at the same frequency of the 
longitudinal counterpart. As a consequence a wider rep¬ 
resentation of the inelastic X-ray scattering data of liq¬ 
uid water was suggested; a Damped Harmonic Oscilla¬ 
tor (DHO) line shape, reflecting the transverse dynam¬ 
ics, was added to the viscoelastic model function which 
accounts for the central peak and the longitudinal dy¬ 
namics [ 13 . A better agreement between experimental 
and fitted line shapes was claimed. 


Quite recently Giordano and Monaco presented new 
experimental results of the spectrum of density fluctua¬ 
tions S{k,uj) of liquid sodium, obtained by inelastic X- 
ray scattering [l^. The novelty and relevance of their 
paper arise from the fact that despite the firm belief that 
liquid metal dynamics was a well settled problem, com¬ 
pletely understood in the generalized Langevin theoreti¬ 
cal framework [ij , they analyze their experimental data in 


terms of a Lorentzian function for the central peak, plus 
two DHO functions for the ^longitudinal and transverse 
excitations". To validate this assumption, they com¬ 
pared the experimental results for the liquid with similar 
findings, derived from the same measurements performed 
on polycrystalline sodium. The correspondence of both 
longitudinal and transverse dispersion curves of liquid 
and polycrystalline phases reinforced the idea that liq¬ 
uid sodium can support transverse acoustic excitations 
beside the well-known longitudinal ones. 

In the present paper we report computer simulation 
results of transverse and longitudinal current correlation 
functions and relative spectra. Our goal is to validate 
(refute) the following experimental outcomes : a) trans¬ 
verse modes, with frequency maxima close to the ones 
obtained by fitting the experimental S{k,ui), are present 
in liquid sodium; b) peak frequencies of these transverse 
modes change with the wavevector k as found experimen¬ 
tally; c) transverse modes become ill defined at increasing 
wavevectors. The definitions of longitudinal and trans¬ 
verse currents and their corresponding correlation func¬ 
tions, CL{k,t) and CT{k,t), are well known Q and there 
is no need to report them here. The Fourier transforms 
CL{k, oj) and Cxik, uj) of the normalized correlation func¬ 
tions will be the subject of our analysis. 

MD simulations have been performed with 500 atoms 
in a cubic box of side L with periodic boundary condi- 
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tions. The well known Price, Singwi, Tosi 1^ potential 
model has been adopted along with the standard Verlet 
algorithm for the integration of the equation of motion 
with a time step dt = 0.256t where St = 0.007ps is the 
unit of time of our system. Two independent runs of 
400.000 steps have been carried out; vectors has been 
chosen of the form k = {n,m,l)/L, performing an av¬ 
erage over all the possible triplets {n,m,l) which real¬ 
ize the same vector amplitude k. The simulations have 
been done at two temperatures, namely T = 395K (just 
above melting) and T = 23K, where the system has been 
brought after a rapid quenching in order to prevent crys¬ 
tallization. The comparison of the results at these dif¬ 
ferent temperatures will give hints on the involved dy¬ 
namical processes. A straight Fourier transform of the 
longitudinal and transverse current correlation functions 
for each k has been applied to obtain the requested spec¬ 
tra. These have been examined by searching the best 
Lorentzian lines fitting the data. The choice of the line 
shape has no particular physical reason in this context, 
but it is a simple tool to trace the maximum of the spec¬ 
trum and its width. We have done the same analysis with 
a Gaussian line shape and no significance differences have 
been observed. 



FIG. 3. Dispersion curve of the longitudinal (A) and trans¬ 
verse (B) current spectra. Triangles are the experimental data 
reported in M- The wavevector k is expressed in A Lines 
are guide for the eyes. 


Typical spectra at the two temperatures and various 
wavevectors are reported in Fig. [T] and Fig. [5] for the 
transverse and longitudinal currents, respectively. A few 


comments are in order to shed light on the different be¬ 
havior at the examined temperatures. First: the different 
value of the transverse spectra at cu = 0 (very close to 
zero at T = 23K, finite at T = 395K). A simple ex¬ 
planation is provided by recognizing that as fc —>■ 0 the 
normalized transverse current correlation function reads: 


CT{k,t) 
Cxik, t = 0) 


exp (-?7 - 1). 

pm 


( 1 ) 


Where rj represent the shear viscosity. As a consequence. 
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( 2 ) 


turns out to be very small at low temperature where the 
shear viscosity becomes extremely large. The same argu¬ 
ments apply at finite k where, in a viscoelastic approxi¬ 
mation, the shear viscosity p is replaced by a wavevector 
dependent viscosity 77 (fc). A more detailed analysis can 
be found in Q and [ig. Second: the appearance, at 
low temperature, of a second peak at larger wavevectors. 
As is evident from Fig. [U when k > 1.13A“^ the trans¬ 
verse spectra show two distinct peaks which are very ev¬ 
ident at T = 23K. It has in fact been shown that, to a 
first approximation, transverse and longitudinal currents 
coincide with the velocity auto-correlation function (see 
Eq. 45 of 0 ) , whose spectrum at low temperature starts 
from zero due to the absence of diffusion processes and 
presents two well clear peaks. This feature is much less 
pronounced in the longitudinal current spectra since it is 
proportional to a fact that suppresses all small fre¬ 
quency characteristics. Third: broadening of the peaks 
as temperature increases. This effect is due to the pres¬ 
ence of more and more transits in the system at higher 
temperature. Such an effect has already been shown in 
the study of density fluctuations and the corresponding 
spectra S{k,u}) [l^. The peak position of the longitudi¬ 
nal and transverse current spectra (dispersion relation) 
are reported in Fig. [2KA) and Fig. |3KB), for the two ex¬ 
amined temperatures along with the experimental find¬ 
ings provided in [l3 |. The transverse dispersion relation 
has been obtained by fitting the spectrum by a single 
Lorentzian line at both temperatures. 

A perfect agreement between experimental and com¬ 
puter simulation is evident for the longitudinal dispersion 
curves. The transverse dispersion curves are more com¬ 
plicated: for 0.2A“^ < k < 1.2A“^ the experimental 
transverse excitations are in reasonable agreement with 
the corresponding simulation spectra. For k > I.2A“^ a 
substantial discrepancy is evident. As anticipated, in or¬ 
der to clarify this behavior we have performed the fitting 
of the spectra by two Lorentzian lines both at low and 
at high temperature. The results are shown in Fig. |4l 
It appears that the experimental results keep track of 
the appearance of two transverse components at higher 
wavevectors. 
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FIG. 4. Same as Fig. |3jB), but fitting done with two 
Lorentzian lines for k > 
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FIG. 5. Wavevector dependence of the current spectra width 
at the two temperatures: (A) Longitudinal, (B) Transverse. 


Finally it is interesting to observe the wavevector de¬ 
pendence of the width of the longitudinal and transverse 
current spectra, when both the spectra are fitted by a sin¬ 
gle Lorentzian line. Results are shown in Fig. [3A) and 
[^B). There is a substantial increase of the width when 
the system is brought to the melting temperature. As 
already discussed in [l^ this result reflects the fact that 
at higher temperature transits occur in larger and larger 
numbers and destroy the coherence of the longitudinal 


and transverse oscillations. Also, at both temperatures 
a linear k dependence of the width is apparent. This is 
an indication that both systems are out of the hydro¬ 
dynamics regime, in which the width would follow a 
dependence. 


In conclusion, MD calculations confirm the presence of 
transverse excitations in liquid sodium for wavevectors 
k > 0.2A“^. Their frequencies are in good agreement 
with those derived from the analysis of the experimen¬ 
tal S{k,u}) reported in [1^. The presence of two peaks 
in the transverse current spectra, as revealed by com¬ 
puter simulation, explain the experimental findings be¬ 
yond I.3A“^. This result is supported by the presence of 
a double peak at low temperature. The large increase of 
the line widths of the longitudinal and transverse exci¬ 
tations at the higher temperature , as derived from MD 
simulations, gives support to the description of the dy¬ 
namical properties in simple liquids in terms of vibrations 
inside random valleys and transits between them. 


We thank G. Monaco for providing the experimental 
data and Pal Jedlovszky for useful comments and sug¬ 
gestions. 
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